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INTRODUCTION 


Models  of  the  magnetic  field  within  the  magnetosphere  have  been 
developed  for  the  purpose  of  numerical  analyses  such  as  particle  tra¬ 
jectory  tracing  and  the  ordering  of  satellite  data.  Roederer  (1969) 
reviewed  some  of  ‘die  early  models.  Within  the  last  few  years,  Olson 
and  Pfitzer  (1974),  Choe  and  Reard  (1974)  and  Mead  and  Fairfield  (1975) 
have  developed  sophist  icated  multiple -term  models  which  include  mag- 
netosphcric  currents  and  the  tilt  of  the  earth's  dipole  axis  with  respect 
to  the  incident  solar  wind  direction.  For  many  purposes,  however,  sim¬ 
ple  models  are  desireable  because  of  their  mathematical  tractability. 

Some  simplified  models  have  been  formulated  from  truncated  versions 
of  the  series  expansions  describing  the  sophisticated  models  (c.  g.  Williams 
and  Mead,  1965,  Olson  and  Pfitzer,  1974). 

In  this  report  an  alternative  simple  model  is  described.  The  basic 
model  is  constructed  from  a  dipole  field  and  a  uniform  field  directed  sun¬ 
ward  in  the  northern  hemisphere  and  antisunward  in  the  southern  hemi¬ 
sphere.  Some  properties  of  the  model  arc  discussed  and  compared  with 
several  other  quantitative  models.  The  incorporation  of  tilt  or  rotation- 
like  features  by  the  addition  of  terms  is  demonstrated  to  illustrate  the 
flexibility  of  the  model. 

THE  BASIC  MODEL 

The  magnetic  field  of  the  present  untilted  magnetospheric  model 


is  described  by  the  equation 


(1) 


B  =  pV  (cos  0/r^)  +  B„p  sin  0  cos  <p  x 


Bx  >  0 

0  <  0  <  it /2 

bt  <  0 

v  /  2  <  0  <  n 

in  which  0  is  colatitude,  0  is  magnetic  longitude  measured  from  mid¬ 
night,  r  is  the  geocentric  distance  and  x  is  the  unit  vector  directed 
perpendicular  to  the  dipole  axis.  The  first  term  in  (1)  describes  the 
dipole  field,  and  the  second  term  describes  a  cur  rent- sheet  field  which 
points  in  opposite  directions  in  the  northern  and  southern  hemispheres. 
The  current  sheet  extends  throughout  the  magnetosphere,  effectively 
replacing  the  magnetopause  and  ring  currents  in  the  inner  magneto¬ 
sphere.  Two  constants,  p  and  B^,,  describe  the  dipole  field  and  the 
current  sheet  field,  respectively.  At  large  distances  from  the  dipole, 
the  field  in  the  magnetotail  is  approximately  equal  to  B^, . 

The  components  of  the  magnetic  field  given  by  eq.  (1)  are: 

3 

Bf  =  -2p  cos  9/r  +  B^,  sin  0  cos  0  (2) 

3 

B0  =  -p  sin  0/r  +  B^,  cos  0  cos  0 

B0  =  -bt  sin  * 

where  B^  changes  sign  at  the  equatorial  plane  as  specified  above. 
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PROPERTIES  OF  THE  BASIC  MODEL 


The  divergence  of  the  field  given  by  eq.  (1)  vanishes,  as  it 
must.  It  was  found  that  the  value  B^,  =  .00015  produces  a  mag¬ 
netosphere  which  resembles  other  quantitative  models  when  the  dipole 
field  strength  p  =  .  31  is  invoked  and  when  r  is  measured  in  earth 
radii.  Several  projected  views  of  the  3-dimensional  model  obtained 
with  these  parameters  are  shown  in  Fig.  1. 

Figure  2  shows  the  field  lines  in  the  noon-midnight  plane  only. 
The  dimpled  magnetopause  of  this  model  is  located  at  r  w  10.  The  co¬ 
latitude  0c  and  radial  disUmce  r  of  the  noon  cusp,  where  the  field 
lines  separate  to  form  the  tail,  can  be  found  by  setting  B  (eq.  (2)) 
equal  to  zero.  The  cusp  occurs  at  the  constant  colatitude 

©c  =  tan*1  sTZ  =  54.7°  (3) 

For  the  parameters  p  =  .  31  and  B.j.,  =  .00015  selected  above,  the 

h  -  0  point  is  located  til  the  radial  distance 

rc  =  [p  tan  0o/Bx]1/3  =  1 3  (4) 

The  magnitude  of  B  along  field  lines  originating  at  various 
dipole  latitudes  is  shown  in  Fig.  3  together  with  similar  data  given 
by  Mead  (1964)  for  his  13-term  model.  Although  the  minimum  B 
ooints  in  the  present  model  lie  enuatorward  rf  fhe  rr>tn?.m~.  in  *hc  Mead 
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The  basic  magnetosphe ric  model  described  by 
eq.  (I).  Projected  views  from  above  the  equa¬ 
torial  plane  (top),  in  the  noon-midnight  plane 
(center),  and  from  the  sun  (bottom)  are  shown. 
Field  lines  are  drawn  for  colatitudes  within  50 
of  the  poles  at  intervals  of  5°  and  for  longitude 
intervals  of  18° 


FRACTIONAL  DISTANCE  ALONG  FIELD  LINE 

-  MEAD  (1964) 

- PRESENT  MODEL 


Fig.  3.  Comparison  of  field  strengths  along  field 
lines  originating  at  various  latitudes  in 
the  present  model  (eq.  (1))  and  the  13- 
term  Mead  model 
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model,  the  behavior  of  B  along  a  field  line  is  similar,  A  qualitative 
comparison  with  several  highly  sophisticated  models  based  on  space¬ 
craft  magnetometer  data  is  illustrated  in  Fig.  4,  which  shows  selected 
field  lines  in  the  noon-midnight  plane.  Aside  from  the  dimple  in  the 
field  lines  near  the  dayside  magnetopause,  the  present  simple  model 
appears  to  agree  with  the  sophisticated  models  about  as  well  as  the 
sophisticated  models  agree  with  each  other  (e.  g.  see  Walker,  1976). 
However,  the  diurnal  variation  at  constant  r  that  occurs  in  these  other 
models  and  in  the  observed  field  does  not  occur  in  the  model  described 
by  eq.  (1 ). 

Selected  field  lines  of  simplified  models  that  are  employed  in 
particle  trajectory  calculations  for  cosmic  ray  access  studies  (Gall 
and  Orozco,  1974)  are  illustrated  in  Fig.  5.  One  obvious  difference 
between  these  models  and  the  present  model  is  the  configuration  of  the 
magnetotail.  In  the  present  model  the  magnetotail  is  composed  of  quasi¬ 
parallel  field  lines  that  gradually  converge  toward  the  equatorial  plane 
where  they  reverse  direction.  Many  of  the  other  simple  models  are 
actually  intended  for  use  only  at  r  <  13  earth  radii  because  they  have 
tail  field  lines  that  diverge  or  show  other  unrealistic  behavior  (e.  g. 
see  Olson  and  Pfitzer,  1974).  Since  the  present  model  more  accurately 
represents  the  observed  magnctotail  (e.  g.  see  Willis  and  Pratt 
1972),  it  is  superior  for  simulating  effects  related  to  the  tail  configuration. 
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Fig.  4.  Comparison  of  field  lines  originating  at  various 
magnetic  latitudes  A  in  the  Olson  and  Pfitzer 
(1974)  (OP  74),  Mead  and  Fairfield  (1975) 

(MF  75),  and  present  magnetospheric  field 
models  in  the  noon-midnight  plane 


A  primary  shortcoming  of  the  present  model  is  the  discontinuity 
of  B  in  the  equatorial  plane.  However,  this  discontinuity  is  easily 
smoothed  in  numerical  analyses  by  introducing  a  factor  tanh  (r  cos  0/ 
6),  where  j  is  a  constant,  in  the  tail  field  term: 

B  =  p  v  (cos  0/r^)  +  B^,  sin  0  cos  0  tanh  (r  cos  0/d)  x  (5) 

(M.  Schulz,  personal  communication).  The  effect  of  this  modification, 
as  shown  by  Fig.  6,  is  essentially  that  of  the  addition  of  a  current  sheet 
of  "thickness"  6  in  the  equatorial  plane.  The  value  of  6  determines  the 
colatitude  of  the  last  "closed"  field  lines  and  the  shape  of  the  magneto¬ 
pause  near  the  equator. 

TILTED  DIPOLE 

The  dipole  in  the  present.  model  can  be  tilted  with  respect  to  the 
tail  by  rotating  the  coordinate  system  for  the  current-sheet  field: 

3 

Br  -  -2p  cos  0/r  -i  B^.  (sin  0  cos  0  cos  a  -  cos  0  sin  a) 

3 

Bq  =  -p  sin  0/r  +  B^,  (cos  0  cos  0  cos  a  +  sin  0  sin  a) 

B0  =  *Bt  Sin  *  cos  a 

(6) 

B^,  >  0  cos  0  >  -  sin  0  cos  0  tan  a 

B^.  <  0  cos  0  <  -  sin  0  cos  0  tan  a 
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where  a  is  the  angle  between  the  dipole  axis  and  the  tail  field  direction 
x.  Selected  field  lines  for  the  tilted  model  with  a  =  30*^  are  shown  in 
Fig.  7.  As  in  other  models  with  a  tilted  dipole  (  e.  g.  ,  Mead  and  Fair- 
field,  1975),  some  field  lines  appear  to  penetrate  the  magnetopause. 
These  lines  could  be  eliminated  if  the  appropriate  magnetopause  cur¬ 
rents  were  included. 

ADDITION  OF  ROTATION 

A  rotation-like  effect  can  be  incorporated  in  the  basic  model  (2) 
by  adding  a  term  to  which  describes  a  cylindrical  vortex  oriented 
along  the  dipole  axis: 


=  - B^,  sin  <f>  -K/r  sin.  0  (7) 

Here  K  characterizes  the  vortex  strength,  or  rotational  velocity. 
Although  this  is  not  a  particularly  realistic  representation  of  the  conse¬ 
quences  of  the  rotation  of  the  magnetosphere,  it  is  qualitatively  con¬ 
sistent  with  the  behavior  that  may  be  found  at  latitudes  above  the  region 
of  corotation. 

Figure  8  illustrates  that  the  effect  of  the  above  addition  to  B^  is 
an  asymmetric  distortion  of  the  field  lines  and  the  magnetopause.  In 
fact,  a  local -time  asymmetry  in  the  magnetopause  shape,  similar  to  the 
distortion  produced  by  the  rotation  term  in  eq.  (8),  has  been  observed 
(Fairfield  and  Mead,  1975,  McDiarmid,  et.  al.  ,  1976). 
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Fig.  7. 


The  titled  ti 

with  >  30° 

shown  in  I* ' i 


Fig.  8.  Illustration  of  the  effect  of  adding  the 
"rotation"  term  -K/r  sin  9  to  "B. 
Equatorial  plane  projections  of  the 
field  lines  used  in  Fig.  1  are  shown 
for  two  values  of  K 


CONCLUDING  REMARKS 


It  is  noteworthy  in  the  present  context  that  Mead  and  Fairfield 
(1975),  in  their  report  describing  a  17-term  model  obtained  by  fitting 
spacecraft  magnetometer  data,  found  that  the  predominant  non-dipole 
term  in  the  series  expansion  of  B  described  a  field  directed  sunward 
in  the  northern  hemisphere  and  antisunward  in  the  southern  hemisphere. 

The  advantages  of  using  the  simple  model  that  has  been  described 
here  are  obvious.  The  step-wise  calculation  of  particle  trajectories  in 
a  model  of  the  magnetospheric  magnetic  field  with  a  realistic  magneto¬ 
tail  can  be  carried  out  using  minimal  computer  time.  Moreover,  the 
effects  on  the  magnetosphere  of  the  spatial  nonuniformity  and  temporal 
variability  of  the  solar  wind  are  simply  modeled  by  adjusting  B^.  For 
example,  asymmetric  distortions  can  be  described  by  B^,  =  B^,  (x,  y,  z). 

The  possibilities  afforded  by  the  parameterization  of  B^,,  together  with 
the  tilt  and  rotation  options,  make  this  an  unusually  flexible  magneto- 
spheric  model  for  both  qualitative  and  quantitative  analyses. 
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LABORATORY  OPERATIONS 


The  Laboratory  Operations  of  The  Aerospace  Corporation  is  conducting 
experimental  and  theoretical  investigations  necessary  for  the  evaluation  and 
application  of  scientific  advances  to  new  military  concepts  and  systems.  Ver¬ 
satility  and  flexibility  have  been  developed  to  a  high  degree  by  the  laboratory 
personnel  in  dealing  with  the  many  problems  encountered  in  the  nation's  rapidly 
developing  space  and  missile  systems.  Expertise  in  the  latest  scientific  devel¬ 
opments  is  vital  to  the  accomplishment  of  tasks  related  to  these  problems.  The 
laboratories  that  contribute  to  this  research  are: 

Aerophysics  Laboratory:  Launch  and  reentry  aerodynamics,  heat  trans¬ 
fer,  reentry  physics,  chemical  kinetics,  structural  mechanics,  flight  dynamics, 
atmospheric  pollution,  and  high- power  gas  lasers. 

Chemistry  and  Physics  Laboratory:  Atmospheric  reactions  and  atmos- 
pheric  optics,  chemical  reactions  in  polluted  atmospheres,  chemical  reactions 
of  excited  species  in  rocket  plumes,  chemical  thermodynamics,  plasma  and 
laser-induced  reactions,  laser  chemistry,  propulsion  chemistry,  space  vacuum 
and  radiation  effects  on  materials,  lubrication  and  surface  phenomena,  photo¬ 
sensitive  materials  and  sensors,  high  precision  laser  ranging,  and  the  appli¬ 
cation  of  physics  and  chemistry  to  problems  of  law  enforcement  and  biomedicine. 

Electronics  Research  Laboratory:  Electromagnetic  theory,  devices,  and 
propagation  phenomena,  including  plasma  electromagnetics;  quantum  electronics, 
lasers,  and  electro-optics;  communication  sciences,  applied  electronics,  semi¬ 
conducting,  superconducting,  and  crystal  device  physics,  optical  and  acoustical 
imaging;  atmospheric  pollution;  millimeter  wave  and  far-infrared  technology. 

Materials  Sciences  Laboratory:  Development  of  new  materials;  metal 
matrix  composites  and  new  forms  ot  carbon;  test  and  evaluation  of  graphite 
and  ceramics  in  reentry;  spacecraft  materials  and  electronic  components  in 
nuclear  weapons  environment;  application  of  fracture  mechanics  to  stress  cor¬ 
rosion  and  fatigue -induced  fractures  in  structural  metals. 

Space  Sciences  Laboratory:  Atmospheric  and  ionospheric  physics,  radia¬ 
tion  from  tlie  atmosphere,  density  and  composition  of  the  atmosphere,  aurorae 
and  airglow;  magnetospheric  physics,  cosmic  rays,  generation  and  propagation 
of  plasma  waves  in  the  magnetosphere;  solar  physics,  studies  of  solar  magnetic 
fields;  space  astronomy,  x-ray  astronomy;  the  effects  of  nuclear  explosions, 
magnetic  storms,  and  solar  activity  on  the  earth's  atmosphere,  ionosphere,  and 
magnetosphere;  the  effects  of  optical,  electromagnetic,  and  particulate  radia¬ 
tions  in  space  on  space  systems. 
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